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PREFACE 


This final report presents results of the research effort supported by the 

Contract NAS8-28249 for the period from June I, 1972 through August 31, 1973 

consisting of three sections listed as follows: 

SECTION IA Non-Linear Pressure Wave Propagation in Pipes 

SECTION IB Attenuation Characteristics of Non-Linear Pressure 
Waves Propagating in Pipes 

SECTION II Experimental Investigation of the Unsteady Forces 
Imposed on a Submerged Orcular Cylinder in 
Accelerating Flow 

Under the Principal Investigator's technical direction, the effort involved in Sections 
IA and B were contributed mainly by C. C. Shih; Section II by David Morrow. To 
make the success of this study possible, valuable assistance was provided by the fol- 
lowing staff members of the Fluid Dynamics Laboratory at the University of Alabama 
in Huntsville: Ronald Scott Dean Montgomery Steve Hecht 

Manfred Loh Ramesh Sanghvi Ming Shih Tom McMeans Barbara Oedamer 

David Walker Cindy Conway Steve K. Mac ey Margaret Cuth ill 

% 

The Instrumentation Laboratory and the Machine Shop of the University of Alabama 
in Huntsville also provided useful assistance to the project. The principal investi- 
gator provided general directions to all phases of the research project. 

Cornelius C. Shih 
Principal Investigator 
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SECTION IA 


NON-LINEAR PRESSURE WAVE PROPAGATION IN PIPES 



SECTION IA 


NON-LINEAR PRESSURE WAVE PROPAGATION IN PIPES 


INTRODUCTION 

Highly intense pressure waves in the near sound field are generated by launch 
vehicles. Characteristics of the propagating acoustical energy in the near sound field 
are highly nonlinear. The use of standard linear scaling techniques for determining 
environmental criteria such as overall sound pressures and power spectral densities, 
has resulted in considerable error. But, to the author's knowledge, there are no ade- 
quate scaling procedures available for use in the near sound field. 

In order to establish a foundation of scaling laws for the highly nonlinear 
waves associated with the launch vehicle, the basic knowledge of the relationships 
among the parameters pertinent to the energy dissipation process associated with the 
propagation of nonlinear pressure waves in thermoviscous media is definitely required. 
Review of literature indicates that there are numerous studies in theoretical aspect by 
Burgers (Ref. I), Hopf (Ref. 2), Cole (Ref. 3), Lighthill (Ref. 4), Rodin (Ref. 5), etc., 
on the Burgers' equation which is generally accepted for approximating one-dimen- 
sional ly the above stated phenomena. However, wave propagation in pipes is a 
typical flow problem solvable by the Burgers' equation with a particular set of 
boundary conditions, if a coefficient of energy dissipation is given. There are several 
studies on the wave propagation in pipes by Blackstock (Ref. 6), Nakamura and 
Takeuchi (Ref. 7), Kantola (Ref. 8), Strunk (Ref. 9), Sawley and White (Ref. 10), 

Brown (Ref. II), D'Souza and Oldenburger (Ref. 12), and Schuder and Binder (Ref. 13). 
These studies are either the treatment of small -amplitude waves, shock waves, plane 
waves of finite amplitude, or inviscid fluids. Some of them were verified experimentally, 
but none of them treats the propagation of nonlinear pressure waves of large amplitude 
through measuring velocity distributions throughout the pipe,in order to gain better un- 
derstanding of the energy dissipation mechanism. 

Specifically, the problem of interest for the present study is to experimentally 
investigate the temporal and spacial velocity profiles of fluid flow in a 3-inch open-end 
pipe of various lengths, produced by the propagation of nonlinear pressure waves for 
various diaphragm burst pressures of a pressure wave generator. As a result, temporal 
and spacial characteristics of wave propagation for a parametric set of nonlinear pressure 
waves in the pipe containing air under atmospheric conditions were determined. 

Velocity measurements at five sections along the pipes of up to 210 ft. in 
length were made with hot-film anemometers for five pressure waves produced by a 
piston. The piston was derived with diaphragm burst pressures at 20, 40, 60 , 80 and 
100 psi in the driver chamber of the pressure wave generator. Experimental results 
ore presented in graphical form and compared with theoretical results obtained by the 
method of characteristics with the isentropic assumption. 
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The Oretlcol Consideration 

The physics of the phenomena considered in this study may be expressed by 
simplified equations of momentum, energy, continuity, and state by following limit- 
ing assumptions: (1) the flow is one-dimensional; (2) the gas process is isentropic; 

(3) the fluid is perfect. 


Momentum: 
Continuity: 
Sonic speed: 


IV + v^v = 

at v ax 


1 dp t ay 

p dt ax 
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p ax 
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Proper rearranging of Eqs. (1), (2), and (3) yields the following characteristic 
equations for the method of characteristics: 



(4) 

left-running 

characteristics 

(5) 



( 6 ) 

right-running 

characteristics 

(7) 


where a denotes the sonic speed, x the axial coordinate along pipe, P the static 
pressure, t the time, y the Specific heat ratio, p the fluid density, V the velocity 
component along x. 

A numerical program was written for temporal and spacial solutions of Eqs. 

(4), (5), (6), and (7). Results of numerical calculations of the flow field are presented 
later for qualitative comparison with the experimental results. 
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Description of Test Facility 

Pressure wove generator 

The pressure wave generator generates pressure waves in test pipes of 3 inches 
in diameter by accelerating an aluminum piston over an 11 -inch distance within the 
beginning section of the test pipe. This generator consists of three major parts: Driver 
Section, 6.6 ft. long; Piston Section, 3.2 ft. long; Test Pipe of various lengths up to 
210 ft. 

I , 

Figure 1 shows the driver section and the piston section. The release of com- 
pressed air from the driver section by electro-thermal ly bursting the diaphragm pro- 
vides the energy to rapidly move the piston and generate the pressure wave in the 
test pipe. Stopping of the piston was accomplished by a double-acting damping system 
and a heavy stopper block. 

The test pipe,as shown in Figure 2,consisted basically of a seamless galva- 
nized steel pipe of 3-inch diameter. A gap of approximately 1/32 inch was provided 
between the aluminum piston and the test pipe. Diameters of the piston and the pipe 
where the piston moves through, are 3.07 and 3.1 inches in order to prevent air leaks 
backward against the piston movement. Five instrumentation ports were installed on 
the test pipe; they were spaced at intervals of 8 and 40 ft. on the upper side of the 
test pipe, respectively for 42 ft. and 210 ft.-pipes. 

A diaphragm bursting technique was developed for relatively low burst pres- 
sures. The diaphragm was made of mylar sheets. The steel wire attached to the 
diaphragm was electrically heated and bursted the diaphragm immediately. This 
method has proved to be reliable for the burst pressures ranging from 20 to 100 psi. 

Instrumentation system 

Five Hot-Film Sensors (Model 1212-20, Thermo-Systems, Inc.) with Constant 
Temperature Anemometers (Model 1033A) were used to measure velocities of the 
unsteady flow of air in the test pipe. The details of sensor installation are shown in 
Figure 4. The relative frequency response of the sensor selected is 40,000 cps; the 
nominal resistance 4-8 ohms. The sensors were Installed at 8 and 40 ft. apart; the 
first sensor was located at 3 and 10 ft. from the pipe entrance; the fifth sensor was 
located at 5 and 40 ft. upstream of the pipe exit; respectively for 42 and 210 ft.- 
pipes. 

Five channels of voltage outputs from the anemometer system were recorded 
on an oscillograph (C£C Datagraph 5-133) at a writing speed of 160 in/sec. The 
sensor calibration curves were established for all sensors, using a sensor calibrator 
(Thermo-Systems Model 1125) along with a precision manometer and a thermometer. 
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Test procedure 


During and after the calibration/ the anemometer outputs were connected 
without an/ alteration to the oscillograph set at the desired recording speed. The 
driver section with a mylar diaphragm was then pressurized to a desired pressure (20, 
40, 60/ 80 or 100 psi), while the piston section was carefully checked for its smooth 
movement within the test pipe and its initial position. As the exact pressure measured 
by a precision pressure gage was reached in the driver section, the diaphragm was 
bursted by electrically heating a circled wire attached to the diaphragm. Approxi- 
mately two seconds before the diaphragm bursting, operation of the oscillograph was 
started. Recording of the voltage outputs from the anemometers was then accomplished. 

Presentation and Discussion of Experimental Data 

The experimental data collected from a series of tests for five burst pressures, 
20, 40, 60, 80 and 100 psi, are presented in graphical form. Sufficient number of 
repetitions, say more than three, was made on each test to ensure the repeatability 
and consistency of the data. 

Figure 3 shows Sensor Positions A, B, C, D, and E across a pipe section. 

Figure 4 presents a typical oscillograph recording of voltage outputs from the five 
sensors located at five sections along the test pipe. 

Temporal velocity profiles at Sensor Position A in five pipe sections for burst 
pressures at 100 and 40 psi for 42 ft, pipe are plotted in Figure 6. Comparison of the 
two sets of profiles indicates that decay of the pressure wave is more notable for the 
case of 40 psi after Section 3. The case of 100 psi seems to Indicate that the wave 
propagation in this length of pipe yields less decay. The acceleration of the propa- 
gated wave recorded at Section 5 for both pressures was caused by the expansion 
wave reflected negatively from the pipe exit. With the propagation speed of approxi- 
mately 1100 ft/sec, the pressure wave took about 10 milliseconds to reflect back from 
the pipe exit to Section 5, causing the flow to accelerate. A slight irregularity of 
the profile at Section 3 has been determined to be strictly instrumental by compara- 
tive tests with other sensors installed at the same section. 

The steep bends occurred on the velocity profile at Section 5 is considered to 
be primarily caused by the wave reflection at the pipe exit. In order to ascertain 
this, an identical test was performed for a 21 -ft. pipe of the same diameter. Test 
results of velocity profiles from the 21-ft. pipe are presented in Figure 7 for com- 
parison with data from the 42-ft. pipe under the same test conditions. The velocity 
profile at Section 3 which is closest to the exit of the 21-ft. pipe is shown to be 
characteristically similar to the velocity profile of Section 5 near the exit of the 
42-ft. pipe. This comparison of data from two different pipe lengths evidently 
illustrates the cause of peculiarity of the velocity profile at the section near the 
pipe exit. 
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Figures 7 , 8, and 9 depict temporal velocity profiles for Sensor Position A 
ot each of five sections respectively, for burst pressures of 40, 60, and 80 psT for 210“ 
ft. pipe. The existence of these phenomena was ascertained by repeoted measure- 
ments under the same test conditions. 


Calculated results of temporal velocity profiles using the method of charac- 
teristics, (Eqs. 4-11), with measured boundary conditions from the piston movement 
were compared with the experimental data for burst pressures of 40 and 80 psi. The 
calculations were carried out with the use of UNIVAC 1108 electronic computer. 
The piston movement was recorded and analyzed with a Fastex high speed movie 
camera at the rate of approximately 4500 frames per second. Since the mathematical 
model was established based on isentropic assumptions, the experiment clearly 
demonstrated that the flow process associated with the wave propagation was not 
isentropic, observing a significant effect of frictional resistance against the flow 
which is subsonic throughout the pipe. As the energy dissipation takes place along 
the pipe, the static pressure tends to decrease along the pipe while the velocity 
decreases rather significantly at a fixed position comparing with the isentropic case. 
These phenomena can be illustrated qualitatively by the following equation modified 
from Equation (1) for unsteady one-dimensional flow in a pipe with friction: 


av r .i ^p + 2 Hi 

at p &x p ax 


- V 


ay 

ax 


(B) 


V 

where t* * T t jwf py: shearing stress along the pipe wall . While the convective 

av 7 

acceleration, V— - is experimentally known to be reasonably small in this case, and 

A X 


av 


a minor positive contribution to ——“is made by the increase in negative pressure 

° 

gradient, a predominant negative contribution to the temporal acceleration, -r— is 

at 


mode by the shearing stress term as demonstrated in both Figures 7 and 9. It is note- 
worthy that a greater deviation of the experimental velocity profile from the calculated 
isentropic data was observed for the case of 40 psi burst pressure, agreeing with the 
theoretical consideration based on Equation (8). Also, it is significant to note that 
the wove fronts of both cases propagated nearly isen tropically as evidenced by the close 
agreement shown in both Figures 7 and 9. 


Figures (10) and (11) present a comparison of numerical calculations with and 
without friction at burst pressures of 40 and 80 psi. The calculations considering friction 
were made with a friction coefficient of f = 0.2. The rate of acceleration is found to 
be subject to the influence of friction. This influence is small in the pipe section near 
the piston; but the effect of friction is growingly pronounced in the direction of propa- 
gation of the pressure wave. 

Figure 10 also demonstrates in a different form these peculiar phenomena of 
velocity distributions across Section 1 and 2 for two time frames respectively, at 100 
and 20 psi burst pressures. Having determined through sensor calibrations thot all hot- 
film sensors were functioning with a reasonable accuracy, say less than 15% in error, 
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the Investigators have determined that the peculiar phenomena may be due to the 
energy dissipation caused by frictional resistance along the wall . The effect of friction 
has been considered predominant to the phenomena. However, further investigations 
and data analysis based on the consideration of energy balance and pressure distri- 
bution are required before mating conclusive statements on these phenomena. 

Figure 11 presents propagation speeds of the pressure waves varying along the 
pipe for five burst pressures. The propagation speed was noted to remain constant for 
most cases with the statistical mean of 1 100 ft/sec and the variance of less than 25 
(ft/see) 3 . No correlation among the axial pipe distance, the burst pressure and the 
propagation speed was obviously noted. ' 

Average temperature and pressure of the test medium, air, were 75°F,and 
29.56 in* Hg absolute under the static condition throughout the test. 

This poper is intended for the presentation of experimental data with a limited 
analysis of the data because of its lack of supporting data, such as pressure distributions 
olong and across the pipe. With some simplifying assumptions, further analysis of the 
data through dimensionless parameters may be conducted. However, the author has 
preferred to present the data in original dimensional form and retained the complete 
data analysis for a later publication. 
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Fig. 1 General View of Driver Section (right). Piston Section 
and Test Pipe (left) of Pressure Wave Generator. 




Fig. 2 General View of Test Pipe. 
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Fig. 11 Comparison of Numerical Calculations with and without 

Friction (f=0.2) at Burst Pressure 80 psi and Sensor Position A. 
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SECTION IB 


ATTENUATION CHARACTERISTICS OF NON-LINEAR PRESSURE 
WAVES PROPAGATING IN PIPES 


INTRODUCTION 

Phenomena of non-linear pressure waves of finite amplitude propagating in the near 
sound field have been of considerable interest in recent years because of their effects 
on the environment. There are numerous studies on the phenomena through the Burgers' 
equation based on one-dimensional approximation/ as noted in the authors' recent 
paper [I ]*. A typical application of these studies is the wave propagation in pipes. 
Recent studies on the wave propagation in pipes have been reported by Blackstock [2], 
Nakamura and Takeuchi [3J, Kantola [4] , Strunk [5], Sawley and White [6], Brown [7], 
D'Souza and Aldenburger [8], Schuder and Binder [9], and Uchida [10]. These studies 
are either the treatment of small-amplitude waves, shock #aves, plane waves of finite 
amplitude, in inviscid or incompressible fluid. Some of them were verified experi- 
mentally, but none of them, to the authors' knowledge, treats the propagation of 
non-linear pressure waves of large amplitude in pipes by measuring velocity distributions 
or profiles throughout the pipe, in order to gain better understanding of the energy 
dissipation mechanism. 

As a step forward in this endeavor, the author [I] reported an experimental investi- 
gation of temporal and spacial velocity profiles of fluid flow in a 3-inch openend pipe 
of various lengths, produced by the propagation of non-linear pressure waves at several 
magnitudes. Experimental results were presented in graphical form and compared with 
theoretical results obtained by the method of characteristics with the isentropic 
assumption. 

This paper presents improved results of the continuing study of propagation phe- 
nomena of the pressure waves in the pipes and provides theoretical explanations of 
peculiar configurations of the velocity profiles foi* the unsteady flow caused by the 
wave propagation. Also, the attenuation characteristics of non-linear pressure waves 
propagating in the pipes were experimentally investigated and analyzed through com- 
parison with two theoretical models, one without friction and one with friction. The 
analysis resulted in a functional relationship among the friction coefficient, flow and 
pipe parameters. 


'Numbers in [ ] brackets refer to references. 
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THEORETICAL CONSIDERATION 

a) Method of Characteristics for Wave Attenuation 

Considering the physics of the phenomena involved in unsteady irreversible 
flow of compressive fluid in general yields the following equations of momentum, 
energy and continuity: 

A7k _ « 


p vft . i 

p t + / dXj u 

P & L'7~ a . . 4 . -j- -$— ( A — — \ 

dXjt + f Pt u 


(1) 

( 2 ) 

(3) 


Where Vj denotes the velocity vector, the stress tensor, H the internal energy, 
e t) the rate of strain tensor, T the temperature , $ the density, and -$the thermal 
conductivity. 

Rearranging of Eqs. (1), (2), and (3) based on one-dimensional approximation with 
linearized diffusion yields a simplified expression as follows, for plane pressure waves 
of finite amplitude, according to Lighthill [11]: 

24> % 3*^ 

t - 1 ax “ ° ax 4 
v 

1 0 


at 


aa . ,, 5a ir-.j 


(4) 

(5) 


Where B denotes the diffusivity defined by Lighthill as 

Where ^is the kinematic viscosity, the dilatational viscosity,^ the dynamic 
viscosity, Pj> the Prandtl number, if the specific heat ratio, the fluid 

particle velocity component along X axis which is coaxial with the pipe, {Ji, the speed 
of an isentropic infinitesimal sound wave at local pressure and density. 

The approximation applied to the theoretical model has made it possible to reduce 
the three fundamental equations to two relatively simple equations. However, some 
cautionary remarks arc essential in regard to the it- applicability for the real fluid flows 
of interest in this paper. 
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The simplified model represented by Eqs. (4) and (5) is assumed to represent plane 
pressure waves of finite amplitude with linearized diffusion effects propagating inX- 
direction. The diffusion effects are linearized by neglecting non-linear terms of low 
order involved in the energy diffusion process due to viscosity and heat conduction as j 
well as relaxation. The relaxation is referred to a process of departure from equilibrium 
distribution of energy in the internal degrees of molecular freedom of, the fluid under 
the wave propagation. In this study, the attenuation of pressure waves is considered 
with the provision that the wave frequency is small compared with the reciprocals of the 
relaxation times of the internal degrees of molecular freedom. 

It may be deduced that attenuation characteristics of the wave propagation are 
attributable to linearized terms of diffusion effects, namely the viscosity and heat 
conduction, as theorized by Lighthill for pressure waves of finite amplitude. 

Examination of Eqs. (4) and (5) indicates that the left-hand sides of the equations 
are the exact forms of the equations of sound waves of finite amplitude under isentropic 
assumption and they constitute the classical theory of shock-wave formation due to 
the wave-form-steepening tendency of the convective terms in the equation of motion, 

• while the right-hand sides are an approximation to the diffusion effects by linearization. 

It should be noted that the approximation will yield a reasonably accurate analysis 
of the phenomena as long as the entropy change in the thermal system is sufficiently 
small, say in the order of at most. denotes the angular velocity in radians 

per second. For a specific application to the pipe flow, the right-hand side of Eq. 

(4) representing the diffusion effects as a cause of attenuation may be modified, 
recognizing that frictional resistance developed along the pipe wall is the predominant 
diffusion effect contributing to attenuation characteristics. Since the frictional re- 
sistance due to wall shear has been expressed conventionally by % 7 TP$^ where 
T 9 denotes the wall shear, J> the pipe diameter, the pipe length, and the inertial 


force is / fhe r i 9 ht-hand side of Eq. (4) may be conveniently equated as, 

S>A 9D U 

Where A * ^-denotes the pipe cross-sectional area. 

For steady viscous flov/s, the wall shear along the pipe, 7 T has been well accepted 


as a function of Darcy-Woisbach's friction coefficient and other flow quantities. 
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( 8 ) 


'r Ljvl 
lo = 8 

Where ^denotes the friction coefficient of Darcy-Weisbachs' head loss 

equation ( — ~~~ ) A^the Reynolds number, the relative roughness of 

inside pipe wall. 

Thus, % may be replaced hy- 'f fr' c t' on coefficient in unsteady flow 
is assumed to be of the same nature as that in steady flow. Steeter [12] is noted to 
apply this approach successfully to waterhammer analysis with nonlinear frictional 
resistance. Eq. (4) then becomes 

. 2<x aa 4V_ „ 

TiT* ^ ax + r - I ax 2D 

Adding and substracting Eqs. (9) and (5) respectively, result in the following 


A + ^ +a ^ L] 
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dx i/ /. 
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1 (y-&r 
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> 2D 
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(13) 


’where 


L i/iy/ -C y* 

replaces -^-“-for purpose of reflecting the change of flow direc- 


tion. Eqs. (10) and (11) represent the left-running characteristics (C 1 ") and Eqs. (12) 
and (13) the right-running characteristics (C ). To solve these equations numerically, 
they are written in finite-difference notation as illustrated in the following figure: 

+ t 


v <«-9 
ft< rf-OU 
tl-l 
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v (*+0 

<X(n+>) 


Fljvtt i. 
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ForEqs. (10) and (Tl) along the C + characteristics, 

VRW-V (n-O +£- { (a.RCn)- afn-o ) 

m —j$[y(>l->)\V(*-0 1 + V'ROOl k'ROO/ ()4) 

— ±-{voi-O+ vR(n > + A<n- 0 +aR(njJ (15) 


For Eqs. (12) and (13) along the C~ characteristics, 

WfA )-\/(n + 0--f2rj-(All(n) - a.cn* i)J 

= - V(m tf/ytino] + VR(n)l WWl] 

s Y-(v(n*i)+ VR (n)+ acn + i)' 1 ’ &#(*>) 

and 

AX *= AXj + A X 2 


06 ) 

(17) 

( 18 ) 


With proper values of and input conditions of the pressure waves specified, 

Eqs. (14), (15), (16) and (17) are solved simultaneously with a sufficient number «f <*X 
and respective time steps which form a grid for numerical calcu lotions of VR and dR 


in a continuous fashion at junctions of the C* and C" characteristics. For completing 
the calculations, proper boundary conditions at both ends of the pipe are required. 

Since the numerical scheme is relatively elementary, details of the numerical procedure 
and computer program to be used with UNIVAC 1107 high-speed computer are omitted 
for the sake of brevity. 

b) Approximate Theoretical Model of Velocity Profiles Across the Pipe Section 
Although Eqs. (1), (2) and (3) are general enough to form a theoretical model 
of cross-sectional velocity profiles along the pipe for unsteady flow generated by the 
pressure wave of finite amplitude, but the exact solution of the model is extremely 
difficult to obtain because of the mathematical complexity and nonlinearity. 

Review of literoture has resulted in finding Uchldas' [ 10 ] work on the pulsating 
viscous flow superposed on the steady laminar motion of incompressible fluid in a 
circular pipe. His study is limited to one-dimensional laminar flow of incompressible 
fluid with superposition of pulsating flows expressed in Lourier series. However, for 
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the purpose of qualitative comparison with experimental results of the cross-sectional 
velocity profiles from the present study, his exact solution is presented below. 




(19) 


Where K^and fc s are constants representing the amplitudes of elemental pulsation, 

I" the radial coordinate in the pipe section, the pipe radius, anc l Jo the 

first solution of Bessel function of zero order. 

Analysis of Eq, (19) shows that calculated results of the velocity profiles exhibit 
a configuration with the maximum velocity near the wall. This finding confirms qualita- 
tively the general configuration of measured velocity profiles in this study. 

EXPERIMENTAL DATA PRESENTATION AND COMPARISON WITH THEORETICAL 
MODELS 

Experimental apparatus including instrumentation systems used for collecting the 
data to be compared with the theory developed have been described in details by the 
author [1], [13], [14] previously. However, for the readers' convenience and continuity 
in understanding of the paper, a brief but pertinent information on the experiment is 
.presented here. 7.62 cm open-end pipes of 15.85 m and 64.0 m in length respectively, 
were tested with hot-film and hot-wire anemometers for velocity profiles across the 
pipe section and along the pipe axis, or centerline of the pipe, of unsteady flows. The 
unsteady flows in the pipes were generated by the piston thrusted into the pipe at various 
speeds which were regulated by the pressure wave generator in terms of diaphragm burst 
pressures at 14.05, 28.10, 42.15, 56.20, and 70.26 g/mm^. The piston speeds were 
measured with a laser device and recorded as a function of time. The air temperature 
was measured at each test for determining the air viscosity. Reynolds numbers of the 
pipe flow ranged for each burst pressure are listed as follows: 


Burst Pressure 
•n g/mm^ (psi) 

14.05(20) 
28.10(40) 
42. 15 (60) 
56.20(80) 
70.26(100) 


Reynolds number 

0~ 41,666 
0~ 62,500 
0 83, 333 

0-97,222 

0 - 111,111 
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The preceding varying Reynolds numbecrsignifies that the unsteady flows generated 
by the propagating waves appear to shift the flow regime from laminar to turbulent. 

With this observation in mind, proper selection of the values of friction coefficients 
for the numerical calculations as well as their comparisons with the experimental 
results was made. Also it was concerned that, since the piston was physically displaced 
only 27.94 cm into the piston in all cases, boundary layer flow on the pipe wall 
seemed unable to develop fully into the pipe axis. But no measurement has been made 
to confirm the above concern in this continuing project of study. Each of all experi- 
mental data was obtained as a result of conducting a statistically sufficient number of 
tests to ensure its' respeatability. 

Typical velocity profiles across the sections 1 and 2 which are 3.60 m and 16.62 m 
downstream from the inlet for burst pressures of 20 and 100 psi were measured and pre- 
sented graphically in Fig. 2. Two half profiles at each section for each burst pressure 
are shown for two time frames (10 and 15 m sec for 100 psi; 20 and 40 m sec for 20 psi). 
All of the profiles show a significant rise in velocity near the wall. This phenomena 
is confirmed by the theoretical results represented by Eq. (19) and also plotted in Fig. 2 
for selected parameters . The values of ^ ^ selected at 10 and 5 are found to be 
approximately equivalent to the cases of burst pressures 100 and 20 psi, respectively. 

It must be emphasize^ however, that the comparison of the experimental and theoretical 
velocity profiles across the section is qualitative because the flow treated in the theoret- 
ical model is assumed to be entirely in the laminar flow regime and the fluid is considered 
incompressible. 

Some typical velocity distributions along the pipe axis measured in the experiment 

A 

for burst pressures of 28.10, 42.15, and 56.20 g/mm z are presented in Figs. 3, 4, and 
5. These velocity distributions are time-dependent velocity curves at the center of 
each of five sections (l through 5) distributed along the pipe at distances of 3.60, 

16.62, 29.45, 42.44, 55.40 m respectively, from the pipe inlet. Included also in 
Figs. 3, 4, and 5 are time-dependent plots of the piston speeds, numerically calculated 
velocity distributions at corresponding locations for respective input conditions of the 
pressure waves. One of these curves in each of Figs. 3, 4, and 5 was calculated with 
the assumption that the entire flow regime in the pipe was fully developed into the 
laminar type where the friction coefficient of steady flow, f - and 
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^**1 .6723 x 10"^ rn^/sec were adopted for the sake of qualitative comparison for 
average air temperature of approximately 80°F during the test. Calculations of other 
velocity distributions are based on the assumptions of fully turbulent flow in the pipe with 
the friction coefficients, f- 0.02, 0.08, and 0.15. The isentropic propagation without 
friction was also calculated and included in Figs. 3, 4, and 5. 

Comparison of the experimental and theoretical curves of time-dependent velocity 
distribution at five sections olong the pipe has resulted in the following observations 
rather than a verification of one by the other: 

a) While repeatability of the experimental curves has been assured through the 
possible best effort, numerical errors developed in theoretical curves was 
found to be less than 1% based on a preliminary error analysis [14]. 

b) Good agreement between the theoretical and experimental curves is noted 
for all three pressures near the inlet region, particular Sections 1 , 2, and 3, 
if the friction coefficient is assumed to be of laminar type or f = 0.02. 

c) The agreement seems to degradate as the waves propagate downstream. Dis- 
agreement is most pronounced at Section 5. However, the use of f = 0.08 in 
the theoretical curve seems to restore some of the agreement. 

d) A perfect agreement between the theory and the experiment is reached as far 
as the speed of wave toe or the front of the pressure waves propagating through 
the pipe, i.e. isentropic propagation. 

e) Theoretically, if all assumptions made are valid, the diffusion effects due to 
viscosity are shown to be compromising with the convective effects which 
steepen the wave form, in a rather linear trend for the laminar friction and 
In a non-linear trend for the turbulent frictioryas far as the time-dependent 
velocity distributions are concerned. The non-linear trend is obviously most 
pronounced for the case of f = 0.15 as shown. 

f) Through the comparison, it may be deduced that the diffusion effects con- 
veniently represented by the friction coefficient are a function of space, time 
and wave magnitude and frequency as well as the relative roughness of the 
pipe for the unsteady flow of interest, 

g) The theoretical model with laminar friction seems to achieve reasonable agree- 
ment with the experiment by using the kinematic viscosity of 0.0016723 m^/sec. 
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Fig. 6 presents a comparison of calculated velocity distributions along 7.62 cm 
(3 in.) and 2.54 cm (I in.) pipes of the same length at the burst pressure of 42.15 
g/mm^ (60 psi). It demonstrates theoretically the effect of pipe diameter on the 
attenuation characteristics under the assumption of f = 0.02. Obviously, the 
reduction in pipe diameter contributes to the non-linear trend of diffusion effects 
as effectively as the increase in the value of friction coefficient. 

The experiment with a I in. pipe of the same length was conducted and its re- 
sults are presented in Figures 7, 8, and 9 for comparison with the numerical calcula- 
tion based on isentropic assumptions qnd with the friction coefficient, f = .02 for 
fully turbulent flow at the three burst pressures. In these I in pipe 1 experiments 
the piston velocities were measured to be slightly higher in comparison with the 
corresponding burst pressures of the 3 in, pipe. However, the effects of diffusion 
as well as the pipe diameter were clearly demonstrated in the I in. pipe experiments 
as shown in the figures. It may be deduced at this juncture that the diffusion 
effects expressed in terms of the friction coefficient, f, are functions of space and 
time for the non-linear propagation of pressure waves in pipes. For the wrought iron 
pipes of 3 in. and I in. in diameter, the friction coefficient, f, may vary from .01 
to .08 and the value of f seems to increase as the waves propagate down the pipe. 
CONCLUSIONS AND RECOMMENDATIONS 

The theoretical models with isentropic assumptions and with friction has been 
constructed and calculated with the use of the method of characteristics. The 
numerical calculations, based on the theoretical models, were compared with the 
experimental results and good agreement between the two was achieved. Through 
' the comparison, the diffusion effects were determined to be a function of space and 
time and they seemed to increase as the waves propagate down the pipes. The 
theoretical models presented in appendices A and B are capable of predicting with 
reasonable accuracy the phenomena of wave propagation in a pipe of a given 
diameter if the friction coefficient can be defined as a function of space and time 
along the pipe. 

However, the experiments performed so far can only provide a limited insight 
into the function of space and time for the friction coefficient. Obviously, a 
continuous effort must be provided to obtain a complete data for the analysis of 


o 
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the function of space and time for the function coefficient, in order that these 
findings comparing the diffusion effects can be made useful for the engineering 
solution of non-linear wave propagations. 
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CHAPTER I 




INTRODUCTION 

Rocket technology's rapid development in the decade following 
the launching of the earth's first artificial satellite brought several 
unsteady fluid flow 1 problems to the attention of the aerospace engi- 
neering community. One area which received attention was the pre- 
diction of the response of flexible, slender bodies to loads imposed 
by ground winds. This was not a new problem, but one with which 
civil engineers had grappled for years in the design of smoke stacks, 
suspension bridges, and tall masts. The aerospace engineer found 
this problem to be somewhat more complex than the civil engineer 
had determined it. The civil engineer usually could make allowances 
for the uncertainty of his design loads by increasing the strength of 
the structure. The same procedure, when followed by the aerospace 
engineer, usually resulted in space vehicles with less than the 
desired capability. Alternate solutions were either to design space 
vehicles with limited regard to the problem or to test aeroelastic 


Unsteady fluid flow is defined as a flow in which the freestream 
velocity is time- dependent in an Eulerian reference frame. 


1 



models of the complex structures to obtain structural responses for 
particular space vehicles in question. 

None of these solutions were entirely acceptable to the aero* 

2 

space engineers, as indicated by the conferences (Refs. 1 through 3) 
held on this subject. One such meeting sponsored by the National 
Aeronautics and Space Administration (NASA) in June 1966 focused 
attention on two aspects of this problem. 

These aspects were: (1) the definition of the atmospheric 
environments surrounding the launch complex and (2) the relation- 
ship between the atmospheric environment and the pressures, 
forces, and moments imposed on the space vehicle. These two 
aspects together with a third, the determination of the response of 
the body once the forcing functions acting on the body are known, 
make up the ground winds problem. The second aspect is the one of 
interest to the fluid dynamicist. The flow field surrounding the 
launch complex and the structural loading conditions arising from 
low-level (altitude < 2, 000 feet) atmospheric movements are illus- 
trated in Figure I- 1. 

Both steady and oscillatory forces are induced on the space 
vehicle by the steady-state ground wind shown in profile. The 
steady force which acts in the direction of the steady wind is 


References are 


located at the end of the main text. 
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FIGURE 1-1. FLOW FIELD SURROUNDING AEROSPACE VEHICLE LAUNCH COMPLEX 
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referred to as drag. A force perpendicular to the direction of the 
steady wind is referred to as lift. The oscillatory force, with com- 
ponents both in and perpendicular to the direction of the steady wind, 
results from a flow phenomenon which occurs over a wide range of 
Reynolds number for flow about a bluff body; i. e. p the space vehicle. 
This flow phenomenon is the alternate shedding of vortices from the 
two sides of the space vehicle. This shedding creates an unequal 
pressure distribution over the vehicle and, hence, oscillatory loads 
are induced on the space vehicle as the alternate vortex shedding 
occurs. Flow deflection caused by the umbilical tower, impingement 
of its turbulent wake, or asymmetries on the vehicle can add a com- 
ponent to the oscillatory forces. Additional components may be 
induced by atmospheric turbulence or gusts. 

The NASA- sponsored meeting (Ref. 1) recommended basic 
research into the unsteady flow phenomenon to establish theories for 
better predicting the time- dependent development of the flow field 
surrounding the launch complex. In response to this recommenda- 
tion, the Fluid Dynamics Laboratory of the University of Alabama 
in Huntsville initiated a series of investigations into this complex 
flow phenomenon. The first research efforts were to establish the 
time-dependent flow field around a rigid, right-circular cylinder 
submerged in an unsteady flow. 



Investigation revealed that for the velocities of interest 
(Voo < 100 ft/sec) the Mach number was low where compressibility 
effects are negligibly small and the dynamic similarity between the 
phenomenon with air and water would justify the use of water as the 
medium for the study. After this determination, the initial effort 
consisted of the design, construction, and calibration of an experi- 
mental hydrodynamic flow facility which could be used to achieve the 
desired unsteady flow. 

The facility, shown in Figure 1-2, consisted of a water tunnel 
equipped with a quick-start mechanism and instrumentation systems 
for the determination of pressure, both on the tunnel walls and on 
the model surface. The water tunnel is a vertically erected steel 
tank 21. 13 feet long and 2. 83 feet in diameter. Two walls running 
the entire tank length help to provide two-dimensional flow but, at 
the same time, reduce the useful test section width to 1.92 feet. 
Plexiglass windows set into the walls allow the model ends to be 
built in to reduce undesirable end effects. The flow visualization 
system also made use of these windows. Location of the test sec- 
tion, together with the dimensions of the various items, are illus- 
trated in Figure 1-3. 

Quick start of the flow is achieved by thermally rupturing a 
mylar diaphragm located at the lower tank end. At the test section, 
water flows having acceleration levels of up to 32. 17 ft/sec 2 can be 
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FIGURE 1-2. OVERALL VIEW OF HYDRODYNAMIC FLOW FACILITY 






produced by varying the opening area of the, mylar diaphragm. 
Operating characteristics of the facility, together with complete 
details of its development and operation, can be found in Reference 4. 

The initial research conducted using this facility consisted 
of the measurement of the time-dependent pressure at various loca- 
tions on the surface of a 3. 0-inch diameter cylindrical model. High- 
speed photography was used in the study of the time-dependent flow 
field surrounding the model. A complete discussion of this initial 
research can be found in Schutzenhofer's thesis, Reference 4. 

A. Primary Objective 

The primary objective of this study was an experimental 
investigation of the unsteady forces imposed on a right-circular 
cylinder immersed in a time -dependent flow. To accomplish this 
objective, it was first necessary to design a system capable of 
measuring the unsteady forces. 

The initial phase in the system design consisted of a survey 
of the available literature to determine the nature (i. e., the fre- 
quency and magnitude) of the unsteady forces. The results of this 
survey are presented in the Bibliography. The amount of experi- 
mental research previously conducted on the unsteady forces 
imposed on a cylindrical model immersed in a time- dependent flow 
was found to be somewhat limited. However, the works of 



Schutzenhofe r (Ref* 4), Sarpkaya (Ref. 5), Pearce (Ref. 6), Jones, 
et al (Ref. 7), Bishop and Hassen (Ref. 8), and Keefe (Ref. 9) were 
most helpful in understanding the phenomena of the unsteady forces 
associated with flows around bluff bodies. 

The attempts to develop a mathematical model of the flow 

phenomena created by time-dependent flows about bluff bodies, of 

which the cylinder is the simplest, have been numerous. The 

classical mathematical approach to this problem was developed by 

von Karman (Ref. 10), and perhaps is still the most significant. 

Others who have formulated mathematical models include Sarpkaya 

(Ref. 11), Phillips (Ref. 12), and Gerrard (Ref. 13). However, 

because of the complex mathematics associated with this problem, 

a formulation which yields results acceptable for general engineering 

use has not been made available to the author's knowledge. 

♦ 

One important question remained unanswered after comple- 
tion of the literature survey: To what degree of accuracy can the 

forces be predicted by two-dimensional inviscid theory? To answer 
this question, the author developed a two-dimensional inviscid 
mathematical model which includes the blockage effects. The deriva- 
tion of this model is presented in Appendix A and results obtained 
from the model are discussed in Chapter V. 



B. Specific Accomplishments 


The specific accomplishments of this study are-: 

(1) Experimental measurement of integrated unsteady 
forces acting on a rigid, right- circular cylinder 
submerged in an accelerating flow 

(2) The design of a simple system that can measure 
the unsteady forces 

(3) Development of a digital data reduction system 

(4) Improvement of the flow visualization system 

(5) Analysis of 16 -millimeter movie films of the 
unsteady Dow phenomena around a circular 
cylinder for various time -dependent Dow 


conditions. 



CHAPTER II 


FORCE MEASUREMENT SYSTEM DESIGN 

The limited experimental data regarding unsteady forces on 
bluff bodies in time-dependent flows left unresolved several impor- 
tant questions concerning the nature of the unsteady forces. These 
questions were: 

(1) Does the sudden decrease in cylinder drag coefficient 
which occurs at about a Reynolds number* of 5 X 10 s 
in steady flow occur in time -dependent flows? 

(2) Is the magnitude of the unsteady forces a function of 
flow acceleration level as predicted by inviscid two- 
dimensional theory? 

(3) Is there a spanwise correlation associated with the 
formation of the vortices at the rear of the cylinder? 

These questions led to the decision to design a system capable 
of measuring the total unsteady forces induced on a cylindrical model 
which spanned the entire width (1.92 feet) of the facility test section. 
This system would measure the integrated effects of the unsteady 


*Based on cylinder diameter, R^q - p 0 VD/p 
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forces acting on the total length of the model. These integrated 
forces could then be used to obtain a spatial average of the unsteady 
lift and drag forces as a function of time. 

/ 

A. Balance Analytical Design 

Several methods of achieving the objective of unsteady force 

measurement were investigated and found feasible. These varied 

from a simple system employing stiff beams in conjunction with 

either load cells or strain gages to a highly complex and expensive 

2 

system utilizing an Inertially Compensated Balance (ICB). From 
the standpoints of simplicity and cost, a system employing a two- 
component strain— gaged beam as the measuring element or balance 
is the most advantageous. The principal disadvantage with this 
type of system, however, is the amplification factor associated with 
the response of the physical system to a dynamic excitation. The 
amplification factor is defined as the ratio of system output to forcing 
function input. 

Assumptions were made that for the range of the design 
forces 2 3 , the force measurement system could be designed to have a 
linear response with small displacements. The characteristics of 


2 For a discussion of the ICB, see References 7 and 14. 

3 The design forces were: lift plane ±100 pounds, drag plane +400 
pounds ♦ 



highly complex systems with small displacements and linear 
responses may be approximated for analytical purposes by combina- 
tions of second-order linear differential equations with constant 
coefficients. Each of these differential equations is assumed to 
represent the equation of motion of a single degree-of-freedom sys- 
tem defined by a natural mode of vibration with a frequency, P n . 
Details of this technique can be found in References 15 and 16. The 
employment of this analysis technique allows the system to be 
analyzed as separate combinations of independent single degree-of- 
freedom systems in both the lift and drag planes. The advantage in 
using this type of analysis is that the inertial forces for those modes 
whose natural frequency is well above the forcing frequency are 

assumed to be small and therefore can be neglected. 

4 

The amplification factor and the phase angle for single 
degree-of-freedom systems are functions of two ratios: 

(1) Ratio of the frequency of forcing function to the 
system's undamped natural frequency 

(2) Ratio of the damping present in' the system to the 
critical damping value for the system. 

An examination of Figure II- 1 shows that for low frequency ratios 
(P/P n < 0.2) and low damping ratios (C/C c < 0.05) the amplification 


4 

Phase angle is defined as the angle between the line. of action of the 
forcing function and the system displacement. 



PHASE ANGLE (deg) AMPLIFICATION FACTOR 
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0.0 1.0 2.0 3.0 

FREQUENCY RATIO, P/P n 


FIGURE II-l. AMPLIFICATION FACTOR AND PHASE ANGLE OF FORCED 
SECOND-ORDER SYSTEM AS A FUNCTION OF FREQUENCY 
RATIOS AND DAMPING RATIOS 



factor and the phase angle are approximately 1.0 and 0 degrees, 
respectively. A system with these characteristics would follow the 
forcing function very closely. 

It would appear that the best force measuring system would 
be the one with the highest natural frequency and the least system 
damping. This is not the case for a system employing a strain- 
gaged beam as the sensing element because the strain of the element 
is an inverse function of the system's natural frequency squared. 

For this reason, the problem of design becomes one of obtaining the 
highest system natural frequency while maintaining an acceptable 
strain for the minimum loads which the system is expected to meas- 
ure. This particular design problem is further complicated by the 
fact that because the model is immersed in a time -dependent flow, 
there is an additional mass of fluid referred to as added mass acting 
with the model as though it were attached to the model. This added 
mass in theory is a function of the model diameter, the test section 
size, and the acceleration level of the fluid (see Appendix A). In 
determining the dynamic response to a physical system, this added 
mass must be combined with the mass of the system. This combined 
mass is usually referred to as the virtual mass of the system. It 
has been shown by Jones (Ref. 17) that the first few natural frequen- 
cies of a beam immersed in a fluid can be predicted numerically by 
conventional methods to within 10-percent accuracy if the virtual 



mass is used in the analysis. Jones also noted that the mode shapes 
associated with these frequencies were not affected by either the 
density or viscosity of the fluid. Since the frequencies of the final 
design were to be determined experimentally, this degree of accu- 
racy was considered sufficient for design purposes. 

Because the magnitude of the added mass quantity is a func- 
tion of acceleration level of the fluid, the force measurement system 
was designed considering the maximum acceleration level which the 
facility could provide; i. e., 32. 17 ft/sec 2 . 

The use of a 3 -inch diameter model limited the cross- 
sectional area available for the balance section, thus limiting the 
stiffness of the balance section. This also had a pronounced effect 
on the frequency- st rain relationship of the balance section. The 
method of force determination from the strain gage output which had 
been chosen (see Section D) depends upon having two strain gage 
bridges in each plane located several inches apart. With the area 
limitation, the distance between the gage locations became an impor- 
tant parameter in establishing the magnitude of the natural frequen- 
cies of the system. 

The various parameters outlined so. complicated the problem 
of designing the system that it was necessary to consider the balance 
and model as a combined cantilevered beam. The effects of varia- 
tions in the kinds of materials, the balance section length, and the 



distribution of cross-sectional area along the length of the combined 
cantilevered beam on the system's natural frequencies were evalu- 
ated using an existing digital computer program (Ref. 18). A con- 
jugate beam analysis of the combined cantilevered beam using the 
maximum anticipated load distributed uniformly was performed to 
determine the ranges of bending stiffnesses required to limit the 
beam's free end static deflection to 0. 05 inch. This range of stiff- 
nesses was then used in determining the initial values of the com- 
puterized parametric analyses. 

The results of the parametric analyses indicated that to 
achieve the high frequency desired (<^200 hertz), it would be neces- 
sary to have a rigid support for the measurement system separate 
from the vertical tank* Additionally, these results indicated that a 
built-up beam of stainless steel and aluminum would give the 
desired frequency while keeping the free end deflection within the 
limits established. 

B. Mechanical Design of Balance and Model 

The final design of the balance consisted of a symmetric 
section of stainless steel with a constant cross-sectional area of 
1.71 in 2 for a length of 3 inches. Matching tapered ends, pre- 
vented from twisting by two 1/4- inch- diameter pins and drawn 
together by a 3/ 8- inch-diameter bolt, were used to join the balance 
section to the support structure and to the model. 



The model final design consisted of three constant cross- 
sectional area sections {two of stainless steel, one of aluminum), 
fastened together by four 1/ 4-inch- diameter dowel pins and a mini- 
mum of 1 inch of Number 32NF threads at each joint. Total weight 
of the model assembly was 10.80 pounds. Engineering drawings of 
each component and the model assembly are presented in Appendix B. 

An installation drawing showing the methods of sealing around 
the model is presented in Figure II-2. The free end of the model was 
allowed to "float" inside a 3. 12 -inch-diameter cutout in the plexi- 
glass window. The model extended 0.25 inch into the window with a 
clearance of 0. 12 inch allowed between the end of the model and the 
window. At the fixed end of the model, one wrap of 1 /2-inch-wide 
flexible cellophane tape was used to seal the gap between the end of 
the model and the adapter section fixed to the facility wall. This 
adapte r section extended 0. 4T inch from the wall so that a gap of 
0.03 inch existed between the model and adapter section. 

The effects of this method of sealing were investigated during 
calibration and no measurable difference between the sealed and 
unsealed cases could be determined for the static loading. Also, no 
effects of sealing were measurable when the frequency check waB 


made. 
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C. Strain Gage Selection and Installation 

The final design of the balance section was based upon the 
use of the BLH Electronics Metal Foil Strain Gage, Number DLB- 
PT-12-2A, in a four-arm active Wheatstone Bridge. These gages, 
made of platinum alloy on a phenolic -glass carrier, were found to 
be most satisfactory considering four major factors: 

(1) Gage factor 

(2) Cost 

(3) Suitability for dynamic strain measurement 

(4) Ease of installation 

Installation of 16 of these gages wired in four (two lift, two 
drag) four-arm active Wheatstone Bridge circuits was accomplished 
according to Figures H-3 and II-4 using Micro-Measurements M- 
Bond 610 epoxy adhesive. Waterproofing of the gage installation was 
accomplished using BLH Barrier E Waterproofing# Figure II-5 is 
a photograph of the waterproof balance section after the attachment 
of a waterproof connector to the leads of the four Wheatstone Bridge 
circuits. 

The four-arm active Wheatstone Bridge was selected because 
of its excellent temperature compensation effects and because it pro- 
vided the maximum electrical output for the same strain level. In 
addition, the fully active bridge has greater accuracy and stability 
than other types of bridges (Ref. 18). 


1 




(a) DISTORTED (UNWRAPPED) VIEW, WIRING FOR 
SIDE-BY-SIDE GAGES 



(b) SCHEMATIC CONNECTION DIAGRAMS 


NOTE: WIRING DIAGRAMS FOR LEFT AND RIGHT SURFACES IDENTICAL 

TO THOSE FOR TOP AND BOTTOM SURFACES 


FIGURE II-3. 


SCHEMATIC DIAGRAMS FOR WIRING WHEATSTONE BRIDGES 
FOUR-ARM ACTIVE, TWO-MOMENT SYSTEM 





0.150 

TYP 



FIGURE 1 1-4. LOCATION 



STRAIN GAGES, TYPICAL OF ALL SURFACES OF BALANCE SECTION 
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D. Method of Force Determination 


Three methods of force determination from Wheatstone 
Bridge output were investigated. These ares 

(1) Basic two-moment system 

(2) Direct force determination 

(3) Composite bridge. 

Each of these methods has its advantages; however, the basic two- 
moment method is more suited to high-speed, digital data processing 
A detailed discussion of the relative merits of each system may be 
found in Reference 19. 

In the two-moment method, two strain gage bridges are 
located on a cantilevered beam some distance apart, as shown in 
Figure II-6. 



FIGURE II -6. CANTILEVERED BEAM WITH TWO STRAIN GAGE BRIDGES 




25 


These bridges determine the beam bending moment at their 
relative locations. It is apparent that the normal force, N, is 
given by 

Ml - M2 
N = — 

whe re 

Ml - beam bending moment at location 1 
M2 - beam bending moment at location 2 
X - distance between bridge locations. 

Thus, if calibration constants can be determined which convert the 
electrical output of the Wheatstone Bridges into bending moments 
and the distance between bridges is known, the amount of force 
acting on the beam can be determined. The outlined procedure is 
applicable to both the lift and drag planes when the normal force, N, 
is given the appropriate definition. 

Initially, both static and dynamic calibration of the force 
measuring system was planned; however, the lack of suitable cali- 
brated test equipment for applying known magnitudes of dynamic 
loads at specific frequencies led to the elimination of the dynamic 
portion of the calibration. The results of an experimental frequency 
determination indicate that the design goal of a high first-mode 
natural frequency was achieved. The experimental value for the 
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first natural-mode frequency was found to be between 168 and 172 
hertz. This frequency was independent of the height of the water in 
the facility. 

Static calibration revealed that design goals of system 
linearity and independency of the lift and drag planes were achieved. 
Details of the procedure used to calibrate the force measuring system 
and calibration results are presented in Appendix C. 



CHAPTER III 

FLOW VISUALIZATION SYSTEM 

One of the objectives of this project was the improvement of 
am existing flow visualization system* This existing system used 
bubbles formed when compressed air was released from hypodermic 
needles connected to a plenum chamber* These spherical bubbles 
were photographed by high-speed cameras during the test runs. 

The problem with this system was that no combination of needle ori- 
fice size 'and plenum operating pressures was found which would 
produce a sufficient supply of small bubbles with approximately uni- 
form diameter. High-speed photography and this technique had been 
employed to obtain some motion pictures of the flow phenomena, but 
improvement in both the number and size of the bubbles was desired 
(Ref. 4). 

Investigation of the hydrogen bubble technique revealed 
that a unique application of this procedure might produce the desired 
improvement. The standard hydrogen bubble technique uses elec- 
trolysis of water to produce hydrogen and oxygen bubbles (Ref. 20). 

A thin platinum wire is used as the cathode and the small hydrogen 
bubbles produced on it are swept away in the passing flow where they 
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are assumed to represent fluid particles in the flow downstream. The 
analyses of this standard technique revealed two problems associated 
with it when used in conjunction with the unsteady flow facility. One 
was the flow disturbance which would be created in front of the model 
by the platinum wire and the other, was the cost associated with the 
wire itself. A known fact is that in steady flow the location of flow 
separation on a circular cylinder is highly dependent on the amount 
of free-stream turbulence in the flow. Thus, the very presence of 
the wire upstream of the model might change the thing being studied; 
i. e,, the flow field around the cylinder. 

The cost of the approximately 3-foot length of aluminum wire 
required to span the facility was prohibitive since the wire section 
probably would be destroyed during each test. The estimated cost 
of replacing the platinum wire after each test exceeded the allotted 
funds for this particular aspect of the study. Therefore, if the 
hydrogen bubble technique was to be incorporated into the flow visu- 
alization system, the two above mentioned problems would have to 
be solved. 

Experimentation with various kinds and diameters of wires 
revealed that a large number of small hydrogen bubbles could be 
produced using a 0.0025-inch diameter music wire as the cathode. 
The number and size of the bubbles produced using straight music 
wire was found to be dependent on the level of electrical current and 
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voltage used. The current and dc voltage levels at which the largest 
number of the most uniform appearing bubbles could be produced 
using the straight 0. 0025-inch-diameter music wire was found to be 
4. 8 amperes of current and 280. 0 volts of dc voltage. 

The problem of the turbulence created by the wire as the flow 
passed was solved by moving the music wire below the model. In 
this application, the voltage was applied before flow initiation. As 
the hydrogen bubbles were generated, they rose from the wire until a 
sheet of bubbles extended from the wire to the surface of the water. 
The water was thus seeded with hydrogen bubbles before the flow 
started. When the flow was initiated, the bubbles were trapped in 
the flow where they were assumed to represent fluid particles. This 
eliminated one source of turbulence found in the standard hydrogen 
bubble technique. 

The two items explained in the preceding paragraphs produced 
most of the desired improvement in the flow visualization system. 
Additional improvement was obtained through the use of a second 
bank, of four 500-watt photo lights, and the lining of the test section 
with 4-foot lengths of household aluminum foil held in place with 
plastic tape. These four items produced a substantial improvement 
in the quality of data obtained using the flow visualization system at 

■o 

a very low cost. Figure III- 1 shows a schematic diagram of the 



improved flow visualization system just before the flow is started. 
Results obtained from the use of this system are presented in 
Chapter V. 
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FIGURE III-l. SCHEMATIC DIAGRAM OF IMPROVED FLOW VISUALIZATION SYSTEM 



CHAPTER IV 


DATA REDUCTION SYSTEM 

A secondary objective of this study was to establish a data 
reduction system which would minimize the problems associated with 
the analyses of the large amount of data obtained for each test run. 
The data reduction system also had to provide a means of reproduc- 
ing the test data bo that simultaneous analyses of data from numerous 
testL runs could be accomplished. Prior experience with other force 
measuring systems convinced the author that recording and retaining 
the data in analog signal form on magnetic tape represented the most 
flexible way of storing and reproducing the test data. The fact that 
data is stored electrically rather than graphically permits automatic 
reduction, allows repetitive recreating of test data, and permits 
time-base expansion or contraction for optimum analysis. 

A. System Components 

An investigation of the recording equipment at the UAH 
Research Institute located the instrumentation necessary not only to 
establish an analog recording system, but also to interface with the 
UNIVAC 1108 digital computing system. The major problem 
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associated with this equipment was that the analog-to-digital conver- 
sion module was limited to basically static data since the maximum 
sampling rate was 15 samples per second. However, the use of one 
of the advantages of magnetic tape recording, time-base expansion, 
permitted the circumvention of this problem. After careful study, 
a digital data reduction system was established using the following 
ins tr umentation : 

(1) Honeywell Medium Band Tape Recorder, Model 7610 

(2) Hewlett-Packard Data Acquisition System, Model 20128 

(3) Honeywell Accudata III Wide Band Amplifiers, 

Model D15A-1003 

(4) System Research Corporation DC Power Supply, 

Model 3560 

(5) System Research Corporation Signal Conditioning 
System, Model 3600 

(6) UNIVAC Digital Computer System, Model 1108 

(7) Stromberg-Carlson Plotter, Model 4020. 

Figure IV- 1 shows a schematic diagram of the data flow 
through this system. Appendix D gives a description of each com- 
ponent together with an explanation of the method of time base expan- 
sion and data point-time correlation. 
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FIGURE IV-1. SCHEMATIC DIAGRAM OF DATA REDUCTION SYSTEM 






B. Method of Filtering 


The method of filtering selected makes use of a recent innova- 
tion in digital computation, the Fast Fourier Transform (FFT). The 
FFT is an algorithm which facilitates signal analysis by means of 
digital computers. It is a highly efficient procedure for computing 
the Discrete Fourier Transform (DFT) of a time aeries. The DFT 
has mathematical properties analogous to those of the Fourier inte- 
gral transform and is particularly useful for power spectrum analy- 
sis and filter simulation on digital computers. Reference 21 
provides a bibliography on the FFT while References 22 and 23 pre- 
sent detailed explanations of this computational tool. 

The FFT digital computer program (Ref. 24) used in this 
analysis was strictly for filter simulation. The digital method of 
filtering was selected over electronic filtering for two reasons. One 
was the excellent roll-off characteristics of a FFT simulated filter 
as indicated by Figure IV-2. This figure represents the frequency 
response of a low pass filter with a cut-off frequency of 20 hertz. 
Figures IV-3 and IV-4 present typical examples of unfiltered and 
filtered drag and lift channel outputs. The low pass filter of Figure 
IV-2 was used in each example and was used throughout the entire 
test series except in a specific filter cut-off frequency analysis. The 
second reason was the man-hours associated with digitizing the test 
data. Electronic filtering of the same test data at different frequen- 
cies required separate digitizing for each frequency because the 
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filtering had to be accomplished with the data in analog form. Each 
test run at one frequency required approximately 4 man-hours to digi- 
tize; thus, electronic filtering represented a highly efficient operation 
since digital filtering required only 2 minutes of digital computing 


time. 



CHAPTER V 


EXPERIMENTAL RESULTS 

These experimental results were obtained using the systems 
discussed in the preceding chapters. Measurements of the integrated 
forces imposed on a 3-inch-diameter, right-circular cylinder sub- 
merged in time -dependent flows were accomplished for two area 
ratios*, 1.0 and 1.4. Originally, the test plan called for force 
measurements in fluid flows established by two additional area 
ratios, 4.0 and 8.35. However, efforts to obtain force measure- 
ments in these flows were unsuccessful because of the small magni- 
tude of the fluid accelerations and velocities produced by these area 
ratios. Several attempts were made to obtain force data for these 

area ratios, especially A 4 /A 2 - 4.0, but the force levels obtained 

2 

were below the resolution sensitivity of the force measuring system. 
Force data obtained for the two area ratios are presented in the 


Area ratio, Aj/A 2 , is defined as the ratio of total test section 
cross-sectional area, no model installed, to exit cross-Bectional 
area. 

Resolution sensitivity is defined as the minimum change in the 
measured variable which produces an effective response in the 
instrument. 
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following sections in both force magnitude and coefficient forms. 
Flow parameters used in nondimensionalizing the force data are 
also presented. A comparison of the measured drag coefficient data 
with steady flow drag coefficient data also is shown with selected 
photographs taken from 1 6 -millimeter film used to illustrate the 
complex flow fields surrounding a circular cylinder immersed in 
a time -dependent flow. Finally, the facility- related and model 
configuration- related parameters which might influence the force 
data are discussed. 

A. Flow Parameters 

Fluid accelerations were computed using static pressures 
measured by flush mounted transducers located at different heights 
in the facility wall. The procedure used in computing fluid 
accelerations from differential static pressure measurements is 
explained by Schutzenhofer (Ref. 4). The equation used in this 
procedure is 


A = g 


A P 

Po A H 


whe re 

A - fluid acceleration 
8 


(V-1) 


acceleration of gravity 



AP - pressure difference between adjacent wall- 
mounted transducers 

p Q - fluid density 

AH - height between the particular pressure transducers* 
Before the fluid accelerations were calculated, the pressure 
measurements were filtered using the low pass digital filter shown 
in Figure IV -3. 

Fluid velocities, V, and fluid displacements, S, were com- 
puted by numerically integrating a segmented, third-order curve 
obtained by “curve fitting' 1 the fluid accelerations. Four different 
segments were required because of the necessity of using data from 
different pairs of transducers as the fluid height decreased. The 
fourth segment of the curve was obtained by extrapolating the fluid 
accelerations obtained for the previous segments. For the area 
ratio 1.0 test runs, this extrapolation became necessary after 0.7 
seconds. The fluid acceleration level computed for all area ratio 
1 . 0 tests did not exceed a 0. 5-percent variation from a 1 -g 
acceleration during the first three segments; consequently, it is 
believed that the use of a constant 1-g acceleration for the fourth 
segment does not significantly influence the area ratio 1.0 coeffi- 
cient data. 

The primary nondimensionalizing factor used in this study 
was that developed by Sc.hutzenhofe r (Ref. 4) which includes the 



fluid acceleration as well as the square of the fluid velocity. This 
method of nondimensionalizing was shown to be appropriate for 
both steady and time-dependent flows. Equation V-2 shows this 
method of nondimensionalizing as it applies to force data using the 
drag force as an example; 


Drag Force 

(0. 5 p Q V z + D A p Q ) * L • D 


(V-2) 


where 

Cj - sectional drag coefficient 
p o - fluid density 

A - fluid acceleration 

V - fluid velocity 

D - model diameter 

L. - model length. 

Typical flow parameters computed for area ratio 1.0 tests 
are presented in Figure V-l. The instantaneous Reynolds number, 
R(t), was calculated as a function of time using the flow velocities, 
V, the cylinder diameter, D, the fluid velocity, p Q , and the fluid 
viscosity, p.. 

Figure V-2 presents the fluid acceleration and velocity 


obtained for a typical area ratio 1.4 test. The Quid 'acceleration 
can be seen to start at a 1-g level, decreasing to zero at approxi- 
mately 0. 78 seconds. This flow data was obtained using the same 
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FIGURE V-l. TYPICAL FLOW PARAMETERS FOR AN AREA RATIO 1.0 TEST 
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procedure as the area ratio 1.0 data. The extrapolation segment 
of the fluid acceleration curve was the zero acceleration period. 

The nondimeneional fluid displacement, S/D, and the 
instantaneous Reynolds number calculated for this typical area ratio 
1.4 test are presented in Figure V-3. 

B. Force Data 

The total lift force measured by the force measurement 
system during a typical area ratio 1.0 test is presented in Figure 
V-4 as a function of nondimens ional fluid displacement, S/D. 
Sarpkaya and Garrison (Ref. 25) have shown that the nondimensional 
fluid displacement is a useful parameter for correlating the effects 
of various magnitudes of time-dependent flows. The sign given to 
the lift force is arbitrary. The large negative lift force indicated 
at an S/D of approximately 66 is a result of the fluid free surface 
nearing the model location. 

Figure V-5 presents a typical drag force measurement for 
an area ratio of 1.0. The sudden decrease in drag force corre- 
sponding to the beginning of the critical Reynolds number regime 
is apparent at an S/D of approximately 14. The sudden increase 
^ drag force at an S/D of approximately 66 is also the result of 
the fluid free surface nearing the model. 
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NONDIMENSIONAL FLUID DISPLACEMENT, S/D 


FIGURE V-4. TOTAL LIFT FORCE AS A FUNCTION OF NONDIMENSIONAL 
FLUID DISPLACEMENT, S/D, FOR A TYPICAL AREA 
RATIO 1.0 TEST 
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FIGURE V-5 . TOTAL DRAG FORCE AS A FUNCTION OF NONDIMENSIONAL 
FLUID DISPLACEMENT, S/D, FOR A TYPICAL AREA RATIO 
1.0 TEST 


i 


Figures V-6 and V-7 present sectional lift {C|)and drag 

'i 

(C^) coefficients calculated using the forces presented in Figures 
V-4 and V-5. Figure V-6 shows the coefficient variations with 
time while Figure V-7 shows the coefficient variations with non- 
dimensional fluid displacement, S/D. 

Figure V-8 shows plots of the lift and drag coefficients 
obtained for area ratio 1.0 tests as a function of the instantaneous 
Reynolds number- The data in Figure V-8 begin at an instantaneous 
Reynolds number of 1 X 10 5 because of the transient dynamics of 
the model-balance combination created at the initiation of the test 
run. Figure V-8a shows the randomness of the lift coefficients 
obtained for all test runs of area ratio 1.0 when plotted as a func- 
tion of the instantaneous Reynolds number. 

A plot of the drag coefficients obtained for the accelerating 
flows produced by the area ratio 1,0 tests is shown on Figure V-8b. 
Although some data scatter is shown below the instantaneous 
Reynolds number of 3 X 10 5 , good agreement is shown in the general 
trend and magnitude of drag coefficient. The variation in minimum 
drag coefficient determined for the area ratio 1.0 flows was from 
0. 184 to 0.230 at instantaneous Reynolds numbers of approximately 
5. 5 X 10 5 . This is a variation of 20 percent. 


^Hereafter referred to as lift and drag coefficients. 
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(a) 



(b) 

FIGURE V-6. LIFT AND DRAG COEFFICIENTS AS A FUNCTION OF 
TIME FOR A TYPICAL AREA RATIO 1.0 TEST 
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FIGURE V- 7. LIFT AND DRAG COEFFICIENTS AS A FUNCTION OF NON- 
DIMENSIONAL FLUID DISPLACEMENT, S/D, FOR A 
TYPICAL AREA RATIO 1.0 TEST 
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The integrated lift force imposed on the circular cylinder 
during a typical area ratio 1.4 test is presented in Figure V-9. 

The bias indicated by the data of Figure V-9 is typical of all lift 
data obtained for area ratio 1.4 tests. No explanation of this bias 
is determined. 

Figure V-10 shows the total drag force determined for the 
typical area ratio 1.4 test. The initial portion, S/D ^ 20, of this 
curve is similar to that obtained for an area ratio 1.0 test, although 
delayed somewhat because of the smaller fluid acceleration and, 
consequently, the fluid velocity. The constant velocity assumption 
for the area ratio 1.4 teBts shows up in that the drag force stays 
at a level consistent with a steady velocity in the critical Reynolds 
number regime after an S/D of approximately 36. 

Plots of the lift and drag coefficients obtained for the area 
ratio 1.4 tests are shown in Figure V-ll. The bias in lift force 
and the randomness of the lift data are clearly shown in Figure V-lla 
while the constant velocity or constant Reynolds number effect is 
apparent in the drag data presented in Figure V-llb. 

The effect of the low pass filter cut-off frequency on the 
drag coefficient data obtained for a typical area ratio 1. 0 test is 
shown in Figure V-12. Although the initial portion (time <0.3 
second) of these data does indicate a more oscillatory pattern as 
the filter cut-off frequency is increased, that portion of the data 
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NONDIMENSIONAL FLUID DISPLACEMENT, S/D 


FIGURE V-9. TOTAL LIFT FORCE AS A FUNCTION OF NONDIMENSIONAL 
FLUID DISPLACEMENT, S/D, FOR A TYPICAL AREA 
RATIO 1 .4 TEST 
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FIGURE V-10. TOTAL DRAG FORCE AS A FUNCTION OF NONDIMENSIONAL 
FLUID DISPLACEMENT, S/D, FOR AN AREA RATIO 1.4 
TEST 
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FIGURE V-12 


EFFECTS OF LOW PASS FILTER CUTOFF FREQUENCY ON DRAG 
COEFFICIENT FOR AN AREA 1.0 TEST 



in the critical Reynolds number regime (time =0.6 second) is 
essentially the same for all cut-off frequencies, 10 hertz through 
45 hertz. The magnitude of the minimum drag coefficient has a 
variation of 10 percent over this filter cut-off frequency range. 

The effects of not including the acceleration term in non- 
dimensionalizing the drag force obtained for an area ratio 1.0 test 
is shown by the data of Figure V-13. These data indicate that the 
acceleration term is important only during the initial portion of the 
test (time < 0. 30 second). However, for unsteady flows with a 
low velocity, the inclusion of the acceleration term might provide 
a more meaningful method of nondimensionalizing test data. 

C . Comparison with Steady Flow Data 

A comparison of the drag coefficients measured during this 
investigation of the forces imposed by unsteady flows with those 
measured by other investigators in steady flows is presented on 
Figure V-14. The conventional practice of nondimensionalizing 
using only the flow velocity (0. 5 p Q V 2 ) was used in this comparison. 
These data agree quite well with those measured by Delaney and 
Sorensen (Ref. 26) on a 4-inch-diameter cylinder. In addition, 
the minimum drag coefficients obtained during this study lie within 
the envelope formed by the minimum drag coefficients determined 
by Delaney and Sorensen and other investigators Such as Jones 
(Ref. 7) and Schmidt (Ref. 27). 
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D. Flow Visualization Data 


A series of six photographs taken from a 1 6 -millimeter movie 
film shot at a rate of 400 frames per second during an area ratio 
1.0 test is presented in Figure V-15. The photographs shown were 
taken from a movie film (Ref, 28) of all the various area ratio tests 
previously discussed. The photographs of Figure V-15, in which 
the flow is from left to right as viewed* show the flow fields sur- 
rounding the right-circular cylinder at instantaneous Reynolds 
numbers ranging from 1.30 x 10 5 to 6. 95 X 10 5 . The instantaneous 
Reynolds number shown for each photograph was obtained by cor- 
relating the number of frames from flow initiation with the instan- 
taneous Reynolds numbers calculated for the constant acceleration 
(1 g) flow produced by the area ratio 1. 0 tests. 

Flow fields surrounding the right-circular cylinder in a 
flow produced by an area ratio 8. 35 test are shown in the seven 
photographs of Figure V-16. These photographs were taken from 
a film shot at a rate of 200 frames per second. The interesting 
feature of these flow fields is the "line" of hydrogen bubbles which 
appears in photographs b through d. This "line" appears in all 
movie films taken of tests at this area ratio in approximately the 
same manner. The "line" briefly appears in some movies of area 
ratio 4. 0 tests; however, it is not nearly as pronounced as in the 
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(e) (f) 
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FIGURE V-15. PHOTOGRAPHS OF FLOW FIELDS SURROUNDING A RIGHT- 

CIRCULAR CYLINDER SUBMERGED IN A 1-g ACCELERATING 
FLOW 
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FIGURE V-16. 


PHOTOGRAPHS OF FLOW FIELDS SURROUNDING A RIGHT- 
CIRCULAR CYLINDER DURING A TYPICAL AREA RATIO 
8.35 TEST RUN 
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area ratio 8. 35 tests. No evidence of this "line" is found in the 
area ratio 1. 0 films. No explanation of the "line" has been 
determined. 

E. Discussion of Results 

This investigation has produced drag coefficient data which 
varies similarly as the well-established steady flow drag coefficient 
variation with Reynolds number in the critical Reynolds number 
range. This similarity is shown in Figure V-14 where the results 
of this investigation and selected steady flow results are plotted as a 
function of instantaneous Reynolds number. The coefficient data 
obtained for these tests lie within the envelope of the steady flow data; 
however, the coefficient data presented herein does exhibit a 20- 
percent variation in magnitude. Besides the usual factors which 
cause experimental force data scatter (i. e. , balance accuracy and 
repeatability), the data presented herein are subject to influences 
which have been shown to especially affect steady flow experimental 
results in the critical Reynolds number range. To the author's 
knowledge, no investigation has proven that these influences should 
not affect time -dependent flow force data in an analogous manner. 

A listing of some of the more prominent influences and a brief dis- 
cussion of their possible effects is given below. These influences 


are: 



(1) Freestream turbulence 

(2) Cylinder surface roughness 

(3) Model end configurations 

(4) Flow three-dimensionality. 

Freestream turbulence has been shown by many investigators 
to be an important factor in defining the drag coefficient variation 
with Reynolds number in steady flow. For example, Surry (Ref. 29) 
recently conducted an extensive investigation. Although these tests 
were conducted in the same facility under supposedly identical con- 
ditions, evidence of different levels of freestream turbulence was 
obtained from close scrutiny of the available movie films. These 
different turbulence levels possibly resulted from the slightly dif- 
ferent bottom opening sequences or the exit diameters which occurred 
during these tests. The different bottom opening sequences and exit 
diameters result from the inability to exactly duplicate the rupturing 
of the mylar diaphragm from test to test. Different opening 
sequences and exit diameters also would have produced the varia- 
tions in acceleration levels measured. Because of the randomness 
of the freestream turbulence, its effects on the force data presented 
herein are difficult to determine and are unknown at this time. 

Cylinder surface roughness has been shown by Achenback 
(Ref. 30) to locally affect the flow around a cylinder submerged in 
steady flow. Foreign matter on the surface of the cylinder was 



controlled during this study by cleaning the smooth model (surface 
finish 16 by 10“ 6 inches, peak to peak) before each test with Methyl- 
Ethyl- Ketone (MEK) and a soft cloth. However, this cleaning did not 
remove surface roughness effects from the force data presented. 

The magnitude of these effects is unknown at this time. 

Model end configurations have been shown by Bearman (Ref. 
31) to locally affect the distribution of base pressure coefficients 
along the span of a circular cylinder submerged in steady flow. 
Bearman notes that local disturbances in the base pressure distribu- 
tion resulted from the different span-wise flow levels which were 
determined to be a function of the model end configuration. Identical 
end conditions were used throughout this study so that all data pre- 
sented should include the same end effects unless other flow param- 
eters coupled with the model end configuration to vary the end 
effects. 

Flow three-dimensionality was a subject of much interest 
during the recent 17th European Mechanics Colloquim as reported by 
Mair and Maull (Ref. 32). One of the Colloquim' s major findings was 
that for a stationary model in steady flow, vortex shedding is not 
really two-dimensional, but varies along the span of the model, with 
changes in phase and possibly amplitude of the shed vortices. How- 
ever, the effects of flow three-dimensionality on the data presented 



herein are unknown at this time since the force measurement system 
was designed to measure the integrated forces along the entire model 
length. 

The facility- related factors affecting experimental data can 
he shown to constitute sufficient error sources to cause a significant 
portion of the previously mentioned 20-percent data variation even if 
the effects of wall interference and model blockage on the force data 
are neglected. Because the investigation of model blockage and wall 
interference effects would have required the construction of addi- 
tional specialized force measurement systems, these items were 
not investigated. Therefore, their contribution to the experimental 
data variation is unknown. An estimate of the error sources attri- 
butable to the other facility- related factors is given below. 

The flow velocity used in obtaining the coefficient data shown 
is estimated to be accurate to within ±3 percent. Since the square 
of the flow velocity is used in nondimensionalizing the force data, 
any error in flow velocity is magnified in the coefficient data pre- 
sented in Figure V-14. This magnification means that a variation 
of approximately 6 percent in coefficient data is possible through 
errors in the velocity determination. 

Accuracy of the force measurement system was not suffi- 
ciently high enough so that some errors were not developed in thfe 



coefficient data presented. The small magnitude of the measured 
forces produced errors which were estimated at approximately 3 per- 
cent based on the results of static calibration. This error does not 
multiply the other error caused by velocity determination, but could 
continue to produce an overall error of approximately 9 percent. 

Thus 9 percent out of the 20-percent coefficient variation can be 
attributed to the method of velocity determination and the accuracy 
of the force measurement system. 

The previous discussions of the parameters and error sources 
affecting the coefficient data presented herein do not invalidate its 
correlation with the steady flow drag coefficient, but rather serve 
to possibly explain the scatter in the unsteady flow drag coefficients. 
The general trend of the drag data appears to be very similar to the 
steady flow variation in this regime of the Reynolds number, 2 X 10 5 
through 7.5 X 10 5 . This trend would not change should all error 
sources combine to give the maximum deviation. In this Reynolds 
number range, the overall velocity effects were found to dominate 
the acceleration effects for the model and facility used in this study. 
Based on the data presented herein, it may be deduced that a rela- 
tionship may exist between the drag coefficients for circular cylinders 
in time-dependent and steady flows. This relationship may be 
obtained if the instantaneous Reynolds number is used in analyzing 
time-dependent flows. However, the Reynolds number range 
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investigated and the number of tests performed in this study do not 
provide a sufficient basis to positively establish a relationship 
between the drag coefficients for right- circular cylinders submerged 
in steady and time -dependent flows. 

The previous paragraphs were concerned primarily with data 
obtained for the drag plane in the Reynolds number range of 2 X 10 5 to 
7. 5 X 10 s , This range represented the region in which the output 
from the data reduction system met the author's sampling criterion. 
Below this range, the balance-model combination was producing high- 
frequency dynamic responses to the flow initiation that exceeded the 
force levels associated with the time-dependent flows. These 
dynamic responses had sufficiently damped at approximately 0.3 
second to allow the data reduction system output to meet the sampling 
criterion. This time corresponds to a Reynolds number of approxi- 
mately 2 X 10 5 . Figure V-12 does indicate, however, that some 
small effects of this dynamic response are present throughout each 
individual test. This can be seen by observing the increasing number 
of oscillatory patterns which appear in the force data of Figure V-12 
as the filter cutoff frequency is increased. Because the Reynolds 
number range below 2X 10 s represents the region which inviscid 
two-dimensional theory could be expected to best predict the drag 
forces, no comparison of the theoretical and experimental results 



were made. To accomplish such a comparison effectively, a faster 
rate of data sampling would have to be used and the measurement 
system's dynamic response effects removed from the force data. 

An unsuccessful attempt was made to correlate the force 
data with the photographic data. The primary reason for this fail- 
ure was that during the period when vortex shedding was well estab- 
lished, the balance dynamic response was dominating the force data. 
When the response had sufficiently damped to allow force determina- 
tion, the wake behind the cylinder had become extremely turbulent 
with no definite frequency of vortex shedding recognizable. 

Comparisons of the data presented herein with data obtained 
by other investigators were attempted. Special efforts were made 
to compare these experimental data with data obtained by the investi- 
gators {Refs. 4 through 9) found to be most helpful in understanding 
the complex flow phenomena. Most of the referenced investigations 
were conducted in steady flow with an oscillating model. However, 
in those cases where a stationary model and time -dependent flow 
were used, either the test Reynolds number was too low or the avail- 
able reports contained insufficient information for the calculation of 
those parameters found important in this investigation. Table V-l 
summarizes the attempts at data comparison with the primary reason 
for failure in each case given. 



TABLE V-l. COMPARISON FAILURE SUMMARY 






CHAPTER VI 

CONCLUSIONS AND RECOMMENDATIONS 

A. Conclusions 

An experimental investigation has been conducted in which 
measurements were made of the unsteady forces imposed on a right- 
circular cylinder submerged in time -dependent flows. Conclusions 
reached as a result of this research effort are as follows: 

(1) The drag coefficients for a circular cylinder in 
time-dependent flows appear to be correlated with 
steady flow drag coefficients in the instantaneous 
Reynolds number range of 2 X 10 s to 7. 5 X 10 s ; i. e . , 
the useful Reynolds number range of this investigation. 

(2) Two-dimensional inviscid theory cannot be used to 
accurately predict the unsteady forces since a varia- 
tion with instantaneous Reynolds number is indicated 
by these results. 

(3) A simple and inexpensive force measurement system 
can be constructed to measure the unsteady loads 
imposed by the time -dependent flows. 
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(4) A means of extracting or eliminating the balance- 
model combination dynamic response effects from 
the force data must be implemented if accurate 
force measurements are to be made at low Reynolds 
numbers in impulsively started time -dependent 
flows. 

(5) A modified version of the hydrogen bubble technique 
can be used to obtain movies of the general qualitative 
features of the complex flow patterns surrounding the 
cylinder. 

B. Recommendations 

The following recommendations are made with regard to 
continuing research into the flow fields surrounding bluff bodies: 

(1) The bibliographical literature search be continued 
to maintain awareness of the efforts and results of 
other investigators. 

(2) Dynamic calibration of the balance-model combina- 
tion used in this analysis be accomplished to verify 
the calibration constants. 

(3) Additional testing of the model configuration used in 
this investigation be accomplished to verify the 
trends of the data presented herein. 
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(4) The effects of model end configurations on the unsteady 
forces be evaluated. 

(5) The effects of cylinder surface roughness on minimum 
drag coefficients be evaluated using one-piece models 
of aluminum machined to exact surface roughness. 

(6) Different diameters of cylinders be tested to determine 
if any model scale or blockage effects can be 
delineated. 

(7) Methods of determining the effects of tunnel wall 
interference be investigated. 
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APPENDIX A. 


THEORETICAL HYDRODYNAMIC MODEL 

A two-dimensional inviscial model was developed for use in 
predicting the test section size and model diameter effects on the 
acceleration-induced forces imposed on a rigid right- circular 
cylinder submerged in a time- dependent flow. To determine these 
effects, the development of two variations of the theoretical hydro- 
dynamic model was necessary. One variation represented an 
accelerating infinite-length cylinder in a bounded stationary flow. 
The other variation represented a stationary infinite-length cylinder 
past which a bounded fluid was accelerated. Figure A- 1(a) shows 
the coordinate system used for the moving cylinder case while 
Figure A- 1(b) shows the coordinate system employed in the station- 
ary cylinder case. 

.Stationary Cylinder Case 

To develop the theoretical model, it was necessary to define 
the complex potential, fl. The complex potential can be determined 
by considering an infinite series of equal- strength doublets located 
equidistant apart along the y axis. The complex potential for this 
arrangement of axes and doublets is given by 
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+oo 

»■ ■ £ £ 

k~ -oo 

where 

p, - moment of the doublet 
b - distance between the doublets 
z - complex coordinate, x + i y 
k - a summation index. 

Recalling that the partial fraction expansion 
k=+oo . 

«*'* = Z rrir (A - 2) 

k : -00 

it follows that 


Z - k bi 


<A-1) 




— tL_ 
2 bi 


cot 


TTZ 

bi 



coth r— 


(A-3) 


The symmetric distribution of the doublets assures that the straight 
lines 


y = (2k + 1) | (k = ... -2, -1, 0, 1, 2...) (A- 4) 

are streamlines. Streamlines have properties which allow their 
replacement by fixed walls (at y = ±kb) without the flow being changed 
in any way. The symmetry of the doublets also allows the superposi- 
tion of a uniform translatory flow with the complex potential, 


«2 = Uz 


(A-5) 



on to the flow created by the doublets. The total complex potential 
for the complete flow is given by 


+ flj — U z + 


2b 


,, 1TZ 
coth ~r~ 
b 


(A-6) 


The complex velocity is defined as the derivative of the com- 
plex potential, «, taken with respect to the complex coordinate, z. 
Performing this operation on Equation A-6 yields, for the complex 
velocity, 


dfl'p 

dz 


u --34- 

2b* 
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sinh 2 


TTZ 
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(A-7) 


Now at the stagnation points, is zero, thus from Equation A-7 


sinh 


2 E* _ , .if E — 
b ~ 2b 2 U 


(A- 8) 


Two stagnation points are located on the real axis equidistant from 
the coordinate system origin. For these two points, the y-axis 
coordinate is zero, thus 


l z l 


b 

IT 



= a 


(A-9) 


where a is the radius of the cylinder represented by a single doublet. 
Therefore, for the condition described. 



which yields for the magnitude of the moment of the doublets. 




2b 2 U 

IT 


sinh 2 


jra 

b 


(A-ll) 


Hence, the equation for the complex potential becomes 

£1t = Uz + — — sinh 2 ~ coth — ~ . (A- 12) 

it b b 

Recalling that 


Slrp — 4 1 ijj 


(A-13) 


where 


<t» - velocity potential 

- stream function . 

Equation A- 12 can be solved for the velocity potential by expanding 
coth irz/b through the use of the following: 

i e 


„ n ,, it r e 

coth -r- = coth 1 v_ — 

b b 


{A- 14) 


and 


coth 
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where the coefficients C n ' are functions of I— J , n=3, 4, 5 . . . . 



Substituting Equation A- 15 into Equation A- 12 yields 
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sinh 2 


/ ira \ b cos 8 
\ b / _ a tr 
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+ Y C n ’ e i(2n * 3)6 . (A- 16). 

T 

where 0 is the angle definition in Figure 4-1. Now using De Moivre's 
theorem and taking the velocity potential, <J>, as the real part of the 
expansion of Equation A- 16 gives 


4> = UX + 


bU 


sinh 2 


(V)f 


b cos 0 . a-ir i0 


air 


3b 


oo 

+ ^ C n ' cos (2n-3) 9 
3 * 


(A- 17) 


Now following Robertson , the force on the cylinder can be obtained 
using the unsteady Bernoulli relationship in terms of an axes system 
moving with a velocity, U. This relationship is given by 

f f + 4 +gh + 9^ -T = II(t) (A- 18) 

where 

dp - pressure differential 

p - fluid density 

g - acceleration of gravity 


Robertson, James M. , "Hydrodynamics in Theory and Application", 
Prentice -Hall, Inc., Englewood Cliffs, New Jersey, pp. 156-158, 

1965 
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h - height above some datum 

q - local flow velocity 

<j> - velocity potential in a fixed frame 

V - velocity of a point in a fixed reference frame in 
terms of the moving axes 

II(t) - an independent function of time. 

For a constant density fluid where gravitational effects are negligible 

and the body is not rotating, Equation A- 18 becomes 


E + 9& _ Za (A- 19) 

p 8t p 


since at large distances, the fluid is moving at a steady velocity 
with a pressure P 0 . 

Now, letting P s represent the steady portion of the pressure 

0 u 

and P u the unsteady pressure (i. e. , that associated with — = 
Equation A- 19 becomes 


dU\ 
dt/ * 



(A-20) 
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Substituting Equation A-21 into Equation B-20 yields 


u = P 


!?!*♦;•«-*(?)[ 


b cos 8 
air 


+ — cos 0 + Yj C n ' cos (2n-3) 0^j | . 

The force, F, due to unsteady motion is defined as the 
gral of the unsteady pressure over the surface of the cylinder 


f = ( p u cos e ds = ( 2ir P u 

J s 


cos 0 a d0 


since 


ds = ad0 - differential surface area of a unit length 
of the cylinder 

0 - the angle defined in Figure A- 1. 

Using Equation A-23 in conjunction with Equation A-22 yields 
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Performing the indicated integration on 0, and letting 
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yields for the force 
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since 



cos (2n-3) 0 cos 8 d0 


0, n greater than 2. 


(A-26) 


Thus the force on a unit length of a stationary infinite -length circu- 
lar cylinder is given by 
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p it a 2 


1 + 


(i**i 


sinh 2 



(A-27) 


Equation A-27 can be used to determine the force on the stationary 
cylinder provided the fluid velocity's time rate of change and the 
geometric characteristics of the model and test section are known. 


Moving Cylinder Case 

The complex potential for this variation of the theoretical 
hydrodynamic model can be obtained by the methods explained in the 
stationary cylinder case (Equation A-6) with the superposition of the 
uniform flow being omitted. This yields 



Following the stationary cylinder case procedure, the velocity 
potential is obtained as 


* - .tab* if) co 8 e + Ec n .cos, 2 n-3,e] 

L 3 

(A-29) 

Again the unsteady pressure's contribution to the force on the 
cylinder can be obtained by using 


Pu - P - P s 



(A-30) 


The use of the velocity potential of Equation A-29 in conjunction with 
Equation A-30 yields 


= P 
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00 1 ) 
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Substituting Equation A-31 into Equation A-23 yields the following 
when the indicated integration is performed 


F -- p tt a 2 


/ 1 b 
£ + ?- 


?■) 5 “ h2 if)} fr 


(A-32) 


Examination of Equations A-32 and A-27 shows a difference 

of a it a 2 — in the indicated force levels. This represents the mag- 
r dt 

nitude of the fluid displaced by the cylinder multiplied by the time 
rate of change of fluid velocity. 
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ENGINEERING DRAWINGS OF MODEL,, 
BALANCE, AND MODEL ASSEMBLY 
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FIGURE B-2. ENGINEERING DRAWING OF BALANCE SECTION 
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FIGURE B-6. ENGINEERING DRAWING OF SUPPORT BALANCE RECEPTACLE 
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APPENDIX C 

CALIBRATION PROCEDURE AND RESULTS 

The employment of Wheatstone Bridges in conjunction with the 
two-moment method of force determination in a force measurement 
system requires calibration to accurately determine force levels. 
This calibration is necessary to determine the proper factors to 
convert the force measurement system's output from volts to 
pounds of force at the model's midpoint* The static calibration pro- 
cedure used in determining the force measuring system conversion 
factors is similar to that used in several wind tunnels to calibrate 
balances used for sting— mounted model applications. The particular 
procedure used assumes that only first-order interaction effects will 
be accounted for. First-order interaction effects result from such 
items as gage-factor variations, variations in gage locations, and 
balance cross-sectional area nonuniformity. Nonlinear interaction 
effects (higher order) result from such items as large deflections 
and mechanical interference. 

Equations C-l through C-4 are first-order interaction equa- 
tions used to determine the bending moments at the four gage loca- 
tions (two lift, two drag) shown in Figure C-l. 
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FIGURE C-l. CANTILEVERED BEAM WITH TWO STRAIN GAGE BRIDGES 
(A AND F) 



where 

Mj - bending moment at aft drag station (A), in-lb 
M 2 - bending moment at forward drag station (F), in-lb 
M 3 - bending moment at aft lift station (A), in- lb 
M 4 - bending moment at forward lift station (F), in-lb 

- drag bridge conversion factor, aft drag station, ~~ 

K 2 - drag bridge conversion factor, forward drag station, 

in- lb 
V 
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Kj - lift bridge conversion factor, aft lift station, — — 

K 4 - lift bridge conversion factor, forward lift station, 

in- lb 

V 

Xj - aft drag bridge output, V 
X 2 - forward drag bridge output, V 
X 3 - aft lift bridge output, V 
X 4 - forward lift bridge output, V 

3M, 2iJ , or a interaction coefficient, where 

j , 2 3 , or 4 indicates either lift 

9 M i,2,3, or 4 or drag plane 


Cursory examination of Equations C-l through C-4 indicates an 
iterative procedure would be necessary to determine solutions if the 
conversion factors were known. Fortunately, two terms in each 
equation are repetitive in this particular application and can be set 
to zero. This results from all gage locations being on a single 
sensing element; therefore, whatever affects one gage affects all 
gages, but in different magnitudes. This is different from conven- 
tional balances which usually have separate elements for sensing 
perpendicular applied force systems. Equations C-5 through C-8 
result from this simplification. 


^ + 3 T^ M * 


M 2 - K 2 x, +|^- M, 


(C-5) 


(C-6) 
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m 3 

= k 3 x 3 

8M i 
0M X 

M x 

(C-7) 

m 4 

= k 4 x 4 

+ *Ma- 
+ 3M 2 

m 2 

(C-8) 


An iterative procedure still is necessary to solve these equations; 
however, the magnitude of the iteration has been drastically reduced. 

The procedure for determining the remaining conversion 
factors is as follows: 

A series of known "dead*' weights is placed in ascending and 
descending order on the model at several known locations and the 
resulting system outputs are recorded. Plots of system output^ versus 
applied force can then be constructed for all gage locations as shown 
in Figure C-2. 

"Best fit" curves can be drawn through the points for each 
load station and their slope, in terms of system output per pound of 
applied force, can be determined. Plots using these slopes and the 
known load stations can then be constructed as shown in Figure C-3. 
From these plots the primary conversion factors, K v through K 4 , 
can be obtained since they are the reciprocals of the "best fit" curves 
shown in Figure C-3. The interaction conversion factors also can be 

^ System output magnitudes will be different between the loaded plane 
and the plane normal to it since the normal plane output represents 
only interaction effects. 
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determined from these plots since they are defined as 


3Mi 

9Mj 


K ik x z jk 


where i and j represent perpendicular planes at the same gage loca- 
tion, k. These nondimensional factors are similar to Poisson's 
ratio in that they relate transverse and longitudinal bending moments 
whereas Poisson's ratio relates transverse and longitudinal strains. 

Figures C-4 through C-8 present the results of carrying out 
the above procedure for the designed force measuring system. The 
final numerical results are presented in Table C-l. These results 
were obtained using the complete data reduction system, including 
filtering, as explained in Chapter IV. Fifty readings were taken at 
each "dead weight" load value indicated in Figures C-6 through C-8. 
The average value of these fifty readings was then used in a "least- 
squares" curve fit digital computer program to obtain the slope values 
shown. First-, second- and third-order curves were obtained 
with the first-order curve showing the smallest standard deviation 
in all cases. The same curve-fit program was used to obtain the 
final slope values of the curves shown in Figures C-4 and C-5. 

It should be noted that the conversion factors for the inter- 
actions are zero. Figure C-8 shows typical results obtained for 
all interaction effects in the first phase of the calibration. These 
results indicate the interaction effects are negligible when compared 


to the primary factors. 



UNIT SYSTEM OUTPUT (V/lb * 10'*) UNIT SYSTEM OUTPUT (V/lb x 10~ 5 ) 
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FIGURE C-4. UNIT SYSTEM OUTPUT AS A FUNCTION OF LOAD STATION - 
DRAG PLANE 


UNIT SYSTEM OUTPUT (V/lb x IQ -3 ) UNIT SYSTEM OUTPUT (V/lb x IQ* 3 ) 
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FIGURE C-6 . DRAG PLANE STATIC CALIBRATION RESULTS 
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FIGURE C-6 - Continued 
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Fifix.SE C-6 - Continued 
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FIGURE C-6 - Concluded 
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FIGURE C-7 - Continued 
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FIGURE C-7 - Continued 
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APPENDIX D 


DATA REDUCTION SYSTEM COMPONENT DESCRIPTION 

The following paragraphs give a description of each of the 
components which make up the Data Reduction System. In addition, 
the methods used to obtain time-base expansion and data point time 
correlation are explained. 

The Honeywell Medium Band Tape Recorder Model 7610 is 
a solid-state electronic recorder with fourteen wide band FM record/ 
reproduce channels capable of recording at seven discrete tape speeds 
from 120 in/sec to 3. 75 in/sec. The FM channels used in this 
analysis have a center frequency of 450 kilocycles (30-percent devia- 
tion) and a band width of 200 kilocycles. Tape flutter, skew, and 
time displacement error are maintained at low levels through the use 
of dual vacuum columns and independent servo control of each reel 
motor. Test data were recorded at a tape speed of 120 in/sec and 
reproduced at a tape speed of 3.76 in/sec. This provided a time- 
base expansion of 32. Therefore, a test run of 1 second was made 
to appear as being of 32-second duration. 

The Hewlett-Packard Data Acquisition System Model 2012B 
is a modular, highly accurate system that accepts analog inputs, 
converts them to digital form, and provides a magnetic tape record 
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of the measurement results which can be input into the UNIVAC com- 
puter. The integrating digital voltmeter which is the heart of this 
system was used on the 1-volt scale with 1 -microvolt resolution. 

The accuracy quoted for this voltage scale is 0.01 percent of the 
reading ±30 microvolts (0. 003 percent full scale). Maximum 
sampling rate for this system when used in the magnetic tape record- 
ing mode is given as 15 samples per second. 

Modification to the data acquisition system was necessary to 
provide accurate means of obtaining time correlation of the test 
data. Since the data acquisition system could digitize only one chan- 
nel at a time, several passes (eight, in most cases) through the 
system were necessary for each test run. Time base correlation was 
achieved by recording a steady-state dc voltage level (5.0 volts) 
during each test run. This voltage level, recorded along with the 
test data, was introduced just before test initiation and removed 
immediately after the test completion. This voltage level was fed 
to the data acquisition system together with each data record and 
was used to start and stop the digitizing /recording modules of the 
system. Since each channel was thus digitized and recorded from the 
same starting signal, data points from all channels could be referred 
to a single time base by use of the digital computer. 

To maintain a check of this method, the triggering dc voltage 
also was supplied to a Beckman timer and the time interval used for 
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each data channel recorded. Results indicated that less than 0.001 
second differential existed between the extremes of the data channel 
time bases. This is less than one-half the time between data points. 
It also was determined that the sampling rate for the modified data 
acquisition system was 14.25 samples per second. The sampling 
rate was investigated over several ranges of data field lengths and 
no measurable variation was determined. 

The combination of the time base expansion and the indicated 
sampling rate produced a maximum sampling rate of 456 samples 
per second (At = 0. 00219 second). The requirement initially estab- 
lished by the author for a digitizing system was 10 samples per cycle 
of dynamic data. This exceeds the criterion generally associated 
with digitizing dynamic data, but should not produce any questions 
about this aspect of the digital data reduction system. Following 
this criterion, the maximum frequency which could be analyzed 
in this test data is 45.6 hertz with all higher frequencies to be 
removed, either by electronic or digital filtering. 

The amplifiers used in this analysis were the Honeywell 
Accudata III D-C Wide Band Differential Amplifier. This is a 
chopper- stabilized, all-transistor amplifier designed for use with 
strain gages and thermocouples . It has a frequency response from 
dc to 20 kilocycles, noise and drift less than 0. 007 percent on typical 



gains, a complete variable gain adjustment with gain stability, and 
linearity better than 99. 99 percent. Because of its stable linear, 
low-noise output, it is suitable for driving a visicorder, multi- 
plexer, analog -to -digital converter, or tape recorder. 

The Signal Research Corporation's dc power supplies and 
signal conditioning system also were used. Figure D-l shows the 
four-wire, calibrate both sides, circuit used to power, calibrate, 
monitor, and balance the strain gage bridges. After static calibra- 
tion, this circuit was used to verify the power settings and balance 
condition before each test run. A similar circuit was also used for 
the three-wall mounted pressure transducers. 

The UNIVAC digital computer and the Stromberg-Carlson 
plotter are not considered basic to the data reduction system since 
the software used in this analysis is applicable to numerous other 
computing and plotting systems. 
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FIGURE D-l . SCHEMATIC DIAGRAM OF STRAIN GAGE BRIDGE CIRCUIT 
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